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Abstract 

In the early 1980's an elementary algorithm for computing conformal maps was discovered 
by R. Kiihnau and the first author. The algorithm is fast and accurate, but convergence was 
not known. Given points zo,...,z n in the plane, the algorithm computes an explicit conformal 
map of the unit disk onto a region bounded by a Jordan curve 7 with z 0l . . . , z n £ 7. We prove 
convergence for Jordan regions in the sense of uniformly close boundaries, and give corresponding 
uniform estimates on the closed region and the closed disc for the mapping functions and their 
inverses. Improved estimates are obtained if the data points lie on a C 1 curve or a K— quasicircle. 
The algorithm was discovered as an approximate method for conformal welding, however it can 
also be viewed as a discretization of the Lowner differential equation. 



50. Introduction 



Conformal maps have useful applications to problems in physics, engineering and mathematics, 
but how do you find a conformal map say of the upper half plane H to a complicated region? 
Rather few maps can be given explicitly by hand, so that a computer must be used to find the map 
approximately. One reasonable way to describe a region numerically is to give a large number of 
points on the boundary. One way to say that a computed map defined on M is "close" to a map to 
the region is to require that the boundary of the image be uniformly close to the polygonal curve 
through the data points. Indeed, the only information we may have about the boundary of a region 
are these data points. 



tThe authors are supported in part by NSF grants DMS-0201435 and DMS-0244408. 
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Figure 1. 



In the early 1980's an elementary algorithm was discovered independently by R. Kiihnau [K] 
and the first author. The algorithm is fast and accurate, but convergence was not known. The 
purpose of this paper is to prove convergence in the sense of uniformly close boundaries, and discuss 
related numerical issues. One important aspect of the algorithm that sets it apart from others: in 
many applications both the conformal map and its inverse are required; this algorithm finds both 
simultaneously. 

The algorithm can be viewed as a discretization of the Loewner differential equation, or as an 
approximate solution to a conformal welding problem. The approximation to the conformal map 
is obtained as a composition of conformal maps onto slit halfplanes. Depending on the type of slit 
(hyperbolic geodesic, straight line segment or circular arc) we actually obtain different versions of 
this algorithm. These are described in Section 1. 

We then focus our attention on the "geodesic algorithm" and study its behaviour in different 
situations. The easiest case is discussed in Section 2: If the data points zq, z\, ... are the consecutive 
contact points of a chain of disjoint discs (see Figures 7 and 8 below), then a simple but very 
useful reinterpretation of the algorithm, together with the hyperbolic convexity of discs in simply 
connected domains (J0rgensen's theorem), implies that the curve produced by the algorithm is 
confined to the chain of discs (Theorem 2.2). One consequence is that for any bounded simply 
connected domain Q, the geodesic algorithm can be used to compute a conformal map to a Jordan 
region O c ("c" for computed) so that the Hausdorff distance between dQ and dQ c is as small as 
desired (Theorem 2.4). 
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In Section 3, we describe an extension of the ideas of Section 2 that applies to a variety 
of domains such as smooth domains or quasiconformal discs with small constants, with better 
estimates. For instance, if <9S1 is a C 1 curve, then the geodesic algorithm can be used to compute a 
conformal map to a Jordan region Q c with dft c € C 1 so that the boundaries are uniformly close and 
so that the unit tangent vectors are uniformly close (Theorem 3.10). The heart of the convergence 
proof in these cases is the technical "self-improvement" Lemmas 3.5 and 3.6. In fact, this approach 
constituted our first convergence proof. 

The basic conformal maps and their inverses used in the geodesic algorithm are given in terms 
of linear fractional transformations, squares and square roots. The slit and zipper algorithms use 
elementary maps whose inverses cannot be written in terms of elementary maps. Newton's method, 
however, converges so rapidly that it provides virtually a formula for the inverses. In Section 4 
we discuss how to apply variants of Newton's method by dividing the plane into four regions, and 
prove quadratic convergence in one region. We plan to address the convergence of the slit and 
zipper variants of the algorithm in a forthcoming paper. 

In Sections 5 and 6, we show how estimates on the distance between boundaries of Jordan 
regions gives estimates on the uniform distance between the corresponding conformal maps to ID>, 
and apply these estimates to obtain bounds for the convergence of the conformal maps produced by 
the algorithm. We summarize some of our results as follows: If d£l is contained in a chain of discs 
of radius < e with the data points being the contact points of the discs, or if d£l is a K-quasicircle 
with K close to one and if the data points are consecutive points on dQ of distance comparable 
to e, then the Hausdorff distance between dfl and the boundary of the domain computed by the 
geodesic algorithm, dQ c , is at most e and the conformal maps ip, ip c onto D satisfy 

sup \tp — ip c \ < Ce p , 

where any p < 1/2 works in the disc-chain case, and p is close to 1 if K is close to one. In the case 
of quasicircles, we also have 

supl^" 1 -^- 1 ! <Ce p 

D 

with p close to one. Better estimates are obtained for regions bounded by smoother Jordan curves. 

Section 7 contains a brief discussion of numerical results. The Appendix has a simple self- 
contained proof of J0rgensen's theorem. 

The first author would like to express his deep gratitude to L. Carleson for our exciting inves- 
tigations at Mittag-Leffler Institute 1882-83 which led to the discovery of the zipper algorithms. 
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1. Conformal mapping algorithms 



The Geodesic Algorithm 



The most elementary version of the conformal mapping algorithm is based on the simple map 
f a : H \ 7 — * H where 7 is an arc of a circle from to a 6 EI which is orthogonal to R at 0. 
This map can be realized by a composition of a linear fractional transformation, the square and 
the square root map as illustrated in Figure 2. The orthogonal circle also meets R orthogonally at 
a point b = |a| 2 /Rea and is illustrated by a dashed curve in Figure 2. 
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Figure 2. The basic map f a 



In Figure 2, c = |a| 2 /Ima. Observe that the arc 7 is opened to two adjacent intervals at with a, 
the tip of 7, mapped to 0. The inverse f~ x can be easily found by composing the inverses of these 
elementary maps in the reverse order. 

Now suppose that z , zi, . . . , z n are points in the plane. The basic maps /„ can be used to 
compute a conformal map of EI onto a region f2 c bounded by a Jordan curve which passes through 
the data points as illustrated in Figure 3. 
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The complement in the extended plane of the line segment from zq to z\ can be mapped onto 
EI with the map 

(p x {z) = iJ- — — 
V z - z 

and (fi(zi) = and (fx(zo) = oo. Set £2 = ^1(22) and if2 = fc 2 - Repeating this process, define 

Cfc = <Pk-l fk-2 O ...O ip-i{z k ) 

and 

fk = fc k ■ 

for k = 2, . . . , n. Finally, map a half-disc to EI by letting 

Cn+i = <Pn • • • <fx{zo) G K 

be the image of zq and set 



2 

z 



1 - ^/Cn+l 



The + sign is chosen in the definition of f n +\ if the data points have negative winding number 
(clockwise) around an interior point of d£l, and otherwise the — sign is chosen. Set 

f = f n+1 o f n O . . . O if 2 O f l 

and 

f' 1 = if' 1 Of' 1 o ...o^. 

Then f~ x is a conformal map of H onto a region S7 C such that zj G <9f2 c , j = 0, . . . , n. The 
portion 7^ of <9S7 C between Zj and Zj + \ is the image of the arc of a circle in the upper half plane 
by the analytic map f^ 1 o . . . o fj 1 . In more picturesque language, after applying ipi, we grab the 
ends of the displayed horizontal line segment and pull, splitting apart or unzipping the curve at 0. 
The remaining data points move down until they hit and then each splits into two points, one on 
each side of 0, moving further apart as we continue to pull. 

As an aside, we make a few comments. As mentioned d£l c is piecewise analytic. It is easy 
to see that it is also C 1 since the inverse of the basic map f a in Figure 2 doubles angles at and 
halves angles at ±c. In fact it is also C2 (see Proposition 3.12). If the data points {zj} lie on 
the boundary of a given region dQ, the analyticity of dQ c also allows us in many situations (see 
Proposition 2.5 and Corollary 3.9) to extend if c analytically across dQ c so that the extended map 
is a conformal map of O onto a region with boundary very close to <9D. Note also that f is a 
conformal map of the complement of Q c , C* \ Q c , onto the lower half plane, C\H where C* denotes 
the extended plane. Simply follow the unshaded region in EI in Figure 3. Finally, we remark that it 
is easier to use geodesic arcs in the right-half plane instead of in the upper-half plane when coding 
the algorithm, because most computer languages adopt the convention — § < arg yfz < \. 

The Slit Algorithm 

Given a region Q, then we can select boundary points z , . . . , z n on dil and apply the geodesic 
algorithm. We can view the circular arcs 7 for the basic maps f a as approximating the image of the 
boundary of Q between and a with a circular arc at each stage. We can improve the approximation 
by using straight lines instead of orthogonal arcs. So in the slit algorithm we replace the inverse of 
the maps f a by conformal maps g a : H — ► H \ L where L is a line segment from to a. Explicitly 

g a ( Z )=C{z- P Y{z + l-p) l - p 

where p = arga/7r and C = \a\/p p (l — p) 1_p . 
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Figure 4. The Slit Maps. 



One way to see that g a is a conformal map, is to note that as x traces the real line from — oo to 
+oo, g a (x) traces the boundary of HI \ L and g a (z) ~ Cz for large z and then apply the argument 
principle. Another method would be to construct Re log g a using harmonic measure as in the first 
two pages of [GM]. As in the basic maps of the geodesic algorithm, the line segment from to a is 
opened to two adjacent intervals intervals on R by f a = g" 1 with f a (a) = and / a (oo) = oo. The 
map f a cannot be written in terms of elementary functions, but an effective and rapid numerical 
inverse will be described in section 4. 

We note that as in the geodesic algorithm, the boundary of the region f2 c computed with the 
slit algorithm will be piecewise analytic. However it will not be C 1 . A curve is called C 1 if the arc 
length parameterization has a continuous first derivative. In other words, the direction of the unit 
tangent vector is continuous. Indeed if g a is the map illustrated by Figure 4, and if is another 
such map then gt, o g a forms a curve with angles 2irp and 2n(l — p) on either side of the curve at 
b = gt, (0). Since analytic maps preserve angles, the boundary of the computed region consists of 
analytic arcs with endpoints at the data points, and angles determined by the basic maps. This 
will allow us to accurately compute conformal maps to regions with (a finite number of) "corners" , 
or "bends". 

The Zipper Algorithm 

We can further improve the approximation by replacing the linear slits with arcs of (non- 
orthogonal) circles. In this version we assume there are an even number of boundary points, 
zq, zi, . . . , z 2n +i- The first map is replaced by 



which maps the complement in the extended plane of the circular arc through zq,Zi, z 2 onto HI. 
At each subsequent stage, instead of pulling down one point £fc, we can find a unique circular arc 

7 




through and the (images of) the next two data points C2/C-1 and Q 2 k- By a linear fractional 
transformation t a which preserves M, this arc is mapped to a line segment (assuming the arc is not 
tangent to M at 0. See Figure 5. 




* * * 

p — 1 p 

Figure 5. The Circular Slit Maps. 



The complement of this segment in H can then be mapped to EI as described in the slit 
algorithm, using g^ 1 where d = a/(l — a/b). The composition h a ^ c = g^ 1 o £ a then maps the 
complement of the circular arc in H onto H. Thus at each stage we are giving a "quadratic 
approximation" instead of a linear approximation to the (image of) the boundary. The last map 
tp n+ i is a conformal map of the intersection of a disc with EI where the boundary circular arc passes 
through 0, the image of Z2 n +i and the image of zq by the composition tp n o . . . o ip 1 . See Figure 6. 
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Figure 6. The Zipper Algorithm. 



If the zipper algorithm is used to approximate the boundary of a region with bends or angles 
at some boundary points, then better accuracy is obtained if the bends occur at even numbered 
vertices {z 2n }- 



Conformal Welding 

The discovery of the slit algorithm by the first author came from considering conformal weld- 
ings. (The simpler geodesic algorithm was discovered later.) A decreasing continuous function 
h : [0, +00) — > (—oo,0] with h(0) = is called a conformal welding if there is a conformal map 
/ of HI onto C \ 7 where 7 is a Jordan arc from to 00 such that f(x) = f(h(x)) for x G K. In 
other words, the map / pastes the negative and positive real half-lines together according to the 
prescription h to form a curve. One way to approximate a conformal welding is to prescribe the 
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map h at finitely many points and then construct a conformal mapping of EI which identifies the 
associated intervals. 

A related problem, which the first author considered in joint work with L. Carleson, is: given 
angles a x , a 2 , • • • , a n and < Xi < x 2 < ■ ■ ■ < x n , find points y n < . . . < yi < so that there 
is a Schwarz-Christoffel map / of H onto a region bounded by a polygonal arc tending to oo with 
angles aj,2n — otj at the jth vertex f(xj) = f{yj)- This map welds the intervals and 
[yj+i,yj], j = 1, . . . ,n. Unfortunately, at the time the best Schwarz-Christoffel method was only 
fast enough to do this problem with polygonal curves with up to 20 bends. 

The basic maps g a can be used to compute the conformal maps of weldings. Indeed, suppose 
yi < < xi, let a = xi/(xi — yi), and apply the map g a {z/{x\ — yi)). This map identifies the 
intervals [yi, 0] and [0, x\\, by mapping them to the two "sides" of a line segment Lcl Composing 
maps of this form will give a conformal map if : EI — ► C \ 7 such that <p([xj, Xj + i]) = (p([yj + i, yj]). 
The final intervals are welded together using the map z 2 . The numerical computation of these 
maps is easily fast enough to compose 10 5 basic maps, thereby giving an approximation to almost 
any conformal welding. Conversely, given a Jordan arc 7 connecting to 00, the associated welding 
can be found approximately by using the slit algorithm to approximate the conformal map from EI 
to the complement of 7. 

The idea of closing up such a region using a map of the form (p n +i was suggested by L. Carleson, 
for which we thank him. 

Since we have been asked about this a couple of times, we note that conformal welding can 
also be defined using the conformal maps to the inside and outside of a closed Jordan curve. If / 
is a conformal map of the unit disc D onto a Jordan region Q and if g is a conformal map of C \ B 
onto C \ fi which maps 00 to 00 then / and g extend to be homeomorphisms of B onto O and C \ B 
onto C\Q respectively. Then the map 

h = f' 1 o g : dB — ► dO 

is a homeomorphism of the unit circle and is also called a conformal welding. Again, if we approx- 
imate a homeomorphism h by prescribing it at finitely many points on the circle, then we can use 
the slit algorithm to identify the corresponding intervals. Simply map the disc to the upper-half 
plane so that the first interval / is mapped to M + , the positive reals, and map the complement of 
the disc to the lower half plane so that desired image h(I) is mapped to 1R + . Apply i^fz and now 
proceed to identify the remaining intervals as above. Conformal welding can also be accomplished 
using the geodesic algorithm. We leave the elementary details to the interested reader. 
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From this point of view, the slit or the geodesic algorithms find the conformal welding of a 
curve (approximately). From the point of view of increasing the boundary via a small curve 7^ 
from zj to Zj+i, the algorithms are discrete solutions of Lowner's differential equation. 

§2. Disc-chains 

The geodesic algorithm can be applied to any sequence of data points zo, z\, . . . , z n , unless the 
points are out of order in the sense that a data point Zj belongs to the geodesic from z^-i to z^, 
for some k < j. In this section we will give a simple condition on the data points zo, z\, . . . , z n 
which is sufficient to guarantee that the curve computed by the geodesic algorithm is close to the 
polygon with vertices {zj}. 

Definition 2.1. A disc-chain Do, Di, . . . , D n is a sequence of pairwise disjoint open discs such 
that dDj is tangent to dDj +1 , for j = 0, . . . , n — 1. A closed disc-chain is a disc-chain such that 
dD n is tangent to dD . 

Any closed Jordan polygon P, for example, can be covered by a closed disc-chain with arbitrar- 
ily small radii and centers on P. There are several ways to accomplish this, but one straightforward 
method is the following: Given s > 0, find pairwise disjoint discs {Bj} centered at each vertex, and 
of radius less than e. Then 

P\(jB j =(jL k 

j 

where {L^} are pairwise disjoint closed line segments. Cover each L k with a disc-chain centered 
on Lk tangent to the corresponding Bj at the ends, and radius less than half the distance to any 
other Li, and less than e. 




Figure 7. Disc-chain covering a polygon. 
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Another method for constructing a disc-chain is to use a Whitney decomposition of a simply- 
connected domain. Suppose £1 is a simply connected domain contained in the unit square. The 
square is subdivided into 4 equal squares. Each of these squares is subdivided again into 4 equal 
squares, and the process is repeated. If Q is a square, let 2Q denote the square with the same 
center, and sides twice as long. In the subdivision process, if a square Q satisfies 2Q C U, then no 
further subdivisions are made in Q. Let U n be the union of all squares Q obtained by this process 
with side length at least 2~ n for which 2Q C f2. If zq G fi, let £l n be the component of the interior 
of U n containing zq. Then dtt n is a polygonal Jordan curve. Note that d£l n consists of sides of 
squares Q with length 2~ n . Thus we can form a disc chain by placing a disc of radius 2~ n /2 at 
each vertex of dfl n . The points of tangency are the midpoints of each square with edge length 2~ n 
on dVt n . 

Yet another method for constructing a disc-chain would be to start with a hexagonal grid of 
tangent discs, all of the same size, then select a sequence of these discs which form a disc-chain. 
The boundary circles of a circle packing of a simply connected domain can also be used to form a 
disc-chain. See for example any of the pictures in Stephenson [SK]. 

If Do, D\, . . . , D n is a closed disc-chain, set 

Zj = dDj n dDj +1 , 

for j = 0, . . . , n, where D n+ i = D . 

Theorem 2.2. If D , D 1 , . . . , D n is a closed disc-chain, then the geodesic algorithm applied to the 
data z , z\ , . . . , z n produces a conformal map ip~ x from the upper half plane EI to a region bounded 
by a C 1 and piecewise analytic Jordan curve 7 with 
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Proof. An arc of a circle which is orthogonal to R is a hyperbolic geodesic in the upper half 
plane H. Let 7^ denote the portion of the computed boundary, d£l c , between Zj and Zj+i. Since 
hyperbolic geodesies are preserved by conformal maps, 7^ is a hyperbolic geodesic in 

For this reason, we call the algorithm the "geodesic" algorithm. 
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Using the notation of Figure 2, each map f' 1 is analytic across R \ {±c}, where / a (±c) = 0, 
and Z" 1 is approximated by a square root near ±c. If is another basic map, then f^ 1 is analytic 
and asymptotic to a multiple of z 2 near 0. Thus o f~ l preserves angles at ±c. The geodesic 7,- 
then is an analytic arc which meets 7 3 -_i at Zj with angle 7r. Thus the computed boundary dil is 
C 1 and piecewise analytic. The first arc 70 is a chord of D and hence not tangent to dD . Since 
the angle at z\ between 70 and 71 is 7r, 71 must enter D\, and so by J0rgensen's theorem (see 
Theorem A.l in the appendix) 

7i C D 1 , 

and 71 is not tangent to dD\. By induction 

j = 0,l,...,n. ■ 

Disc-chains can be used to approximate the boundary of an arbitrary simply connected domain. 

Lemma 2.3. Suppose that ft is a bounded simply connected domain. If e > 0, then there is 
a disc-chain D , . . . , D n so that the radius of each Dj is at most e and dfl is contained in an 
e-neighborhood of UDj . 

Proof. We may suppose that Vl is contained in the unit square. Then for n sufficiently large, the 
disc chain constructed using the Whitney squares with side length at least 2~ n , as described above, 
satisfies the conclusions of Lemma 2.3. ■ 

The Hausdorff distance dn in a metric p between two sets A and B is the smallest number 
d such that every point of A is within p-distance d of B, and every point of B is within p-distance 
d of A. The ^-metrics we will consider in this article are the Euclidean and spherical metrics. 

A consequence is the following theorem. 

Theorem 2.4. If Q is a bounded simply connected domain then for any e > 0, the geodesic 
algorithm can be used to find a conformal map f c of D onto a Jordan region Vl c so that 



d H (dn,dn c ) < e, 
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(2.1) 



where dn is the Hausdorff distance in the Euclidean metric. If is a Jordan curve then we can 
find f c so that 

sup|/(z) - f c (z)\ < e, 

z€V> 

where f is a conformal map of D onto f2. 



Proof. The first statement follows immediately from Theorem 2.2 and Lemma 2.3. To prove the 
second statement, note that the boundary of the regions constructed with the Whitney decompo- 
sition converges to dO, in the Frechet sense. By a theorem of Courant [T, page 383], the mapping 
functions can be chosen to be uniformly close. ■ 



We note that if 0, is unbounded, Lemma 2.3 and Theorem 2.4 remain true if we use the 
spherical metric instead of the Euclidean metric to measure the radii of the discs and the distance 

to on. 

There are other ways besides using the Whitney decomposition to approximate the boundary 
of a region by a disc-chain and hence to approximate the mapping function. However, Theorem 
2.4 does not give an explicit estimate of the distance between mapping functions in terms of the 
geometry of the regions. This issue will be explored in greater detail in Sections 5 and 6. 

The von Koch snowflake is an example of a simply connected Jordan domain whose boundary 
has Hausdorff dimension > 1. The standard construction of the von Koch snowflake provides a 
sequence of polygons which approximate it. By Theorem 2.4 the mapping functions constructed 
from these disc-chains converge uniformly to the conformal map to the snowflake. 




Figure 8. Approximating the von Koch snowflake. 
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It is somewhat amusing and perhaps known that a constructive proof of the Riemann mapping 
theorem (without the use of normal families) then follows. Using linear fractional transformations 
and a square root map, we may suppose Q is a bounded simply connected domain. Using the 
disc-chains associated with increasing levels of the Whitney decomposition for instance, Vl can be 
exhausted by an increasing sequence of domains Vl n for which the geodesic algorithm can be used 
to compute the conformal map f n of fl n onto B with f n (zo) = and f' n {zo) > 0. Then by Schwarz's 
lemma 

, x , fm(w) 

u n (w) = log 

fn(w) 

for n = m + 1, m + 2, . . . is an increasing sequence of positive harmonic functions on fi m which 
is bounded above at zq by Schwarz's lemma applied to J" 1 , since ft is bounded. By Harnack's 
estimate u n is bounded on compact subsets of Vl and by the Herglotz integral formula, log j^^jy 
converges uniformly on closed discs contained in Q m . Thus f n converges uniformly on compact 
subsets of fl to an analytic function /. By Hurwitz's theorem / is one-to-one and by Schwarz's 
lemma applied to /~ , / maps Q onto D. 

In the geodesic algorithm, we have viewed the maps f c and (p~ l as conformal maps between EI 
and a region Vl c whose boundary contains the data points. If we are given a region fi, and choose 
data points {z k } G dfl properly, then the next proposition says that the computed maps (p c and 
ip~ x are also conformal maps between the original region fl and a region "close" to H. 

Proposition 2.5. If D , . . . , D n is a closed disc-chain with points of tangency {zk}, and if ft is a 
simply connected domain such that 

n 
k=0 

then the computed map ip c for the data points {z^}® extends to be conformal on fl. 

We remark that changing the sign of the last map <£> n +i m the construction of (p c gives a 
conformal map of the complement of the computed region onto M. We choose the sign so that the 
computed boundary winds once around a given interior point of fl. 

Proof. Without loss of generality Q D Ufc=o an d hence dfl C UdDk- The basic map f a in 
Figure 2 extends by reflection to be a conformal map of C* \ (7 U 7^) onto C* \ [— c, c], where 7^ 
is the reflection of 7 about R. In fact, if a C H is any connected set such that 0, a G a then f a is 
conformal on C* \ (a U a R ), where a R is the reflection of a about R. In particular, if U is a simply 



connected region contained in EI with 0, a £ dU, then C \ f a (U U U R ) consists of two open sets 
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V U — V where (0, c) G V and (0, — c) G — V. Here f/^ denotes the reflection of the set U about R 
and —V = {—z : z G V}. 
Set 

V'fe = ^fc ° • • • ° fi 

and 

W fc =Vfc(C*\{D U...UD n }). 
Then we claim C* \ {W k U W^} consists of 2(n + 1) pairwise disjoint simply connected regions: 

n 2k 

c* \ {w k u W*} = (J MDj) u ^(d*)* u (J u kJ , 

j=k j=l 

where each region U k ,j is symmetric about R and R C Lij^Ukj. The case A; = 1 follows since 
V'i(C* \ Dq) is bounded by two lines from to oo. As noted above, the image of il)k-i{D k -i) U 
ipk-i(Dk-i) R by the map (p k consists of two regions V and —V. The claim now follows by induction. 
Note 5 k = ip k (dQ n dD k -i) is a simple curve connecting a point — c k to and to c k , satisfying 
z G <5fc if and only if — z G 5^, since \J z 2 + 1 is odd. Since each 99^ extends to be one-to-one and 
analytic on Vfc-i(-Cfc-i) and since t/ n ,j, J = 1, ...In are disjoint, the map ip n is one-to-one and 
analytic on By direct inspection, the final map (p n +i extends to be one-to-one and analytic, 
completing the proof of Proposition 2.5 ■ 

As one might surmise from the proof of Proposition 2.5, care must be taken in any numerical 
implementation to assure that the proper branch of V z 2 + c 2 is chosen at each stage in order to 
find the analytic extension of the computed map to all of Q. 

§3. Diamond-chains and Pacmen 

If we have more control than the disc-chain condition on the behavior of the boundary of a 
region, then we show in this section that the geodesic algorithm approximates the boundary with 
better estimates. We will first restrict our attention to domains of the form C \ 7 where 7 is a 
Jordan arc tending to 00. 

Definition 3.1. An e-diamond D(a,b) is an open rhombus with opposite vertices a and b 
and interior angle 2e at a and at b. If a = 00, then an e-diamond D(cc,b) is a sector {z : 
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I arg(z — b) — 9\ < e}. An e-diamond-chain is a pairwise disjoint sequence of e-diamonds 
D(zq, zi), D(z 1 ,Z2), ■ ■ ■ D(z n -i, z n ). A closed e-diamond-chain is an e-diamond-chain with z n = 
zo- 

See Figure 9. Let B(z,R) denote the disc centered at z with radius R. 
Definition 3.2. A pacman is a region of the form 

P = B(z , R)\{z:\ arg(A(z - z ))\ < e}, 
for some radius R < oo, center zq, opening 2e > 0, and rotation A, |A| = 1. 

Let C\ be a constant to be chosen later (see Lemma 3.7), and let zq = oo. 

Definition 3.3. We say that an e-diamond-chain D(oo, z\), D(z 1 ,Z2), ■ ■ ■ D(z n -\, z n ), satisfies the 
e-pacman condition if for each 1 < k < n — 1 the pacman 

P k = B(z k , R k ) \{z:\ arg(— — — J | < e}, 

^Z k — Zk+l' 

with radius R k = C\\z k+ i — z k \/e 2 satisfies 




The pacman P k in Definition 3.3 is chosen to be symmetric about the segment between z k and 
z k+ i with opening 2e equal to the interior angle 2e in the diamond-chain. Note that the e-diamond 
D(z k -i, z k ) may intersect P k . 




Figure 9. A Diamond-chain and a Pacman. 
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When z = oo, the first map in the geodesic algorithm is replaced by <pi(z) = \\fz — ~z[. The 
argument of A can be chosen so that ^1(^2) is purely imaginary, in which case the boundary of the 
constructed region contains the half- line from Z2 through z\ and 00. We will henceforth assume 
that 

D(oo,Zi) = {z : |arg( - — — )| < e}. 

\zi - Z2J 

Theorem 3.4. There exist universal constants e > and Ci such that if an e-diamond-chain 

D(oo, z 1 ),D(z 1 ,z 2 ), D(z n - 1 , z n ) 
satisfies the e-pacman condition with e < £q , and if 

Zk + l — Zk 



arg 



Zk — Zk-l 



< (3-D 



for k = 2, . . . , n — 1, then the boundary curve 7 computed by the geodesic algorithm with the data 
zq = 00, z\ , . . . , z n satisfies 



C (J (D(z k -!,z k ) U{z k }\ 

k=l ^ ' 



7 

Moreover, the argument 6 of the tangent to 7 between z k and z k +i satisfies \9— arg(zk+i — z k )\ < 3e. 

To prove Theorem 3.4, we first give several lemmas. 

Lemma 3.5. There exists Eq > such that if e < £q, and if Q, is a simply connected region bounded 
by a Jordan arc d£l from to 00 with 

{z : I arg z\ < ir — e} C O, 

then the conformal map f of H + = {z : Rez > 0} onto $7 normalized so that /(0) = and 
/(oo) = 00 satisfies 

|argz 2 /'U)|<|, (3.2) 

where z = / _1 (1). 

The circle C Zo which is orthogonal to the imaginary axis at and passes through zq has a 
tangent vector at z with argument equal to 2argz - The quantity argZQ/'(z ) in (3.2) is the 
argument of the tangent vector to f(C Zo ) at f(zo). 
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Proof. We may suppose that \z \ = 1. Set 

/(*) 



g(z) = log • 



z 2 



Then |Im<7(z)| < e on <9IH + and hence also on H + , and | argz | < §, since /(^o) = 1- Set a = ^ 
and 

A = e 09 ^ - 1 = ^ 2q , 

and 

1 + z 

t(z) = Re^o + ilmz - 

Then r is a conformal map of B onto M + such that r(0) = z and 9? is a conformal map of the strip 
{|Imz| < e} onto ID so that tp(g{zo)) = 0. Thus h = ip o g o t is analytic on B, bounded by 1 and 
h(0) = 0, so that by Schwarz's lemma 



\ V '(g(z o ))\\g'(z o )\\r\0)\ = \h'(0)\ < 1. 



Consequently 



and hence 



f(*o) 



/(zo) ^0 
|arg^/'(z )| 



| 9 ^)|<?£g^< * 



arg z + arg 



<-+sin- 1 



7rRezo tt cos | ' 
zof'(zo) 



f(zo) 

£ 



IT COS 



- .(> + I) s + 0(A 



This proves Lemma 3.5 if e is sufficiently small. 



Lemma 3.6. Let f2 satisfy the hypotheses of Lemma 3.5. Ife< Eq/2, then the hyperbolic geodesic 
7 from to 1 for the region Q lies in the kite 

P = {z : I argz| < e} n {z : | arg(l - z)\ < —}, 

and the tangent vectors to 7 have argument less than |e. 



Proof. By J0rgensen's theorem, 7 is contained in the closed disc through 1 and which has slope 
tane at 0. Likewise 7 is contained in the reflection of this disc about R and hence | argz| < e on 
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7- This also shows that 7 is contained in a kite like P but with angles e at both and 1. In the 
proof of Theorem 3.4, however, we need the improvement to ^ of the angle at 1. 

By Lemma 3.5, a portion of 7 near 1 lies in P. Suppose w\ G 7 with | argwi| = 5 < e and 
then apply Lemma 3.5 to the region -^Q with e replaced by e + 5. Then the tangent vector to 7 
at wi has argument where 

5 

|0-argu/i| < -(e + |argu;i|). (3.3) 

Since | aigwi\ < e, we have \9\ < |e. Moreover (3.3) also implies 6 < |e, when arg^i < 0, and 
6 > — |e when argu>i > 0. But if wi is the last point on jddP before reaching 1, this is impossible. 
Thus 7 C P, proving the lemma. ■ 



The next lemma improves Lemma 3.5 by only requiring that the portion of <9f2 in a large disc 
lies inside a small sector. 

Lemma 3.7. There is a constant C\ so that ifs<Eo/2 and if dQ is a Jordan arc such that € dfl, 
dnn{\z\ > C\/e 2 } ^0, and 

c 

{z : I arg z\ < tt — e and \z\ < — } C il, 
then the conformal map f : M+ — ► with /(0) = and |/(oo)| > satisfies 



argzg/'^)! < -, (3.4) 



where zq = f 1 (1). 



Proof. As before, we may assume \z \ = 1. Let w(z, E, V) denote harmonic measure at z for En V 
in V \ E. Set R = and E# = 5(0, i?) = {|z| < i?}. Then by Beurling's projection theorem and 
a direct computation 



uj{1, 0B R , n) < uj{1, 0B R , B R \ [-R, 0]) = - tan" 1 ( -L) . 
By the maximum principle 



(3.5) 



/(*>) 
ar g — — 



<e + (27r + e)-tan _1 ( j 



< 



lie 
10 ' 
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for C\ sufficiently large. Since f(z ) = 1, 

lie 

\^gz \<—. 

Next we show that there is a large half disc contained in fl Br). Set 

S = ini{\w — ilmzol : w £ H + and f(w) £ dBj^}. 

Using the map 

z — ilmzo — S 



z — ilmzo + S 

of H + onto B and Beurling's projection theorem again, 

u(z ,f- 1 (dB R ),m + ) > u(zo,[S,oo) +ilmz ,M + ). 
Then by (3.5), (3.6) and an explicit computation 



-tan _1 f -^y ) > -tan" 1 ! ; 5^ V 



For e sufficiently small, this implies 

i 

5(o,^)nH+cr 1 (fini?(o,§-)). 

Now follow the proof of Lemma 3.5 replacing r with a conformal map of B onto EI + n{|z| < 

such that r(0) = z - Then r'(0) = 2Rez + O(-t-) and for Ci sufficiently large, (3.4) holds. 

c? 



Following the proof of Lemma 3.6 (replacing 5/6 by 9/10), the next corollary obtains. 

Corollary 3.8. Suppose dQ is a Jordan arc such that £ dQ, dO, n {\z\ > Ci/e 2 } / 0, and 

C 

{z : | argz| < tt — £ and \z\ < — — } C Q. 

£ z 

If £ < £q/2, then the hyperbolic geodesic 7 from to 1 for the region lies in the kite 

9e 

P = {z : I argz| < e} n {z : | arg(l — z)\ < —}■ 

Moreover, the tangent vectors to this geodesic have argument at most 3e. 
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Proof of Theorem 3.4. Use the constant C\ from Lemma 3.7 in Definition 3.3. As in Theorem 
2.2, let 7j denote the portion of the computed boundary d£2 between zj and Zj+i- By construction 
7o U 71 is a half line through zq = oo, z\, and zi- Make the inductive hypotheses that 

fc-1 fc-1 

J 7 , C \jD( Zj ,z J+1 ) (3.8) 

3=0 3 = 

and 

7fc _inP fc = 0. (3.9) 

Since the e-diamond chain D(oo, zi), D(z\, z<i), ■ ■ ■ D(z n -i, z n ) satisfies the e-pacman condition, 
(3.8) and (3.9) show that the hypotheses of Corollary 3.8 hold for the curve 7 = Uq _1 7j and hence 
7fc C D(zk, Zfc+i). Also by Corollary 3.8 and (3.1), 

7fc n P k+1 = 0. 

By induction, the theorem follows. ■ 

If the hypotheses of Theorem 3.4 hold, then the proof of Proposition 2.5 gives the following 
Corollary. 

Corollary 3.9. If O and the diamond chain D(zk,Zk+i) satisfy the hypotheses of Theorem 3.4, 
then the conformal map ip c computed in the geodesic algorithm extends to be conformal on Q, U 

\Jk=0 D ( z k,Zk+l)- 

The next Theorem says that for a region 17 bounded by a C 1 curve, the geodesic algorithm 
with data points zo, zi, . . . , z n produces a region O c whose boundary is a C 1 approximation to dQ. 

Theorem 3.10. Suppose Q is a Jordan region bounded by a C 1 curve d£l. Then there exists 
So > 0, depending on dQ so that for 5 < So, dQ, is contained in a closed S -diamond-chain D = 
UD(zk, Zk+i) and so that dfl c , the boundary of the region computed by the geodesic algorithm, is 
contained D. Moreover if ( G dVL c and if a G dVL with \( — a\ < S then 

\v<-V a \<M, (3-10) 
where tj^ and r] a are the unit tangent vectors to <9f2 and <9f2 c at ( and a, respectively. 

Proof. There were two reasons for requiring that zo = 00 in Theorem 3.4. The first reason was to 
assure that 

(Uo"V)n (C\B(z k ,R k ))^fH (3.11) 
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as needed for Lemma 3.7. The second reason is the difficulty in closing the curve, since Lemma 3.7 
does not apply. The difficulty being that a pacman centered at z n will contain zq if z is too close to 
z n . Since d£l G C 1 , we may suppose that the 5-diamond chain D(zo, z\),D{z\, Z2), ■ ■ ■ , D(z n -i, z n ) 
satisfies the pacman condition. Note that this requires z n to be much closer to z n _ 1 than to z . 
Since dfl G C 1 , if |z n — zo| is sufficiently small, we can find two discs 

n-l 

A p c C\ (J D(z k ,z k+1 ), 


for p = 1,2, with 

{z , z n } = 9Ai n dA 2 c Ax n A 2 c L>(z n , z ), 

where D(z n ,zo) is a 5-diamond. By J0rgensen's theorem, as in the proof of Theorem 2.2, the 
geodesic 7„ from z n to zq is contained in Ai n A 2 . Then by the proof of Theorem 3.4, dtt c is 
contained in the 5-diamond chain. The statement about tangent vectors now follows from Corollary 
3.8. ■ 



We say that {z k } are locally evenly spaced if 



D 



< 



Zk — Zk-l 



z k — Zk+1 



(3.12) 



for some constant D < 00. Note that the spacing between points can still grow or decay geometri- 
cally. We define the mesh size /x of the data points {zj} to be 



V-{{zj}) = sup|z fc - z k+ i\. 
k 

We say that a Jordan curve T in the extended plane C* is a i^T-quasicircle if for some linear 
fractional transformation r 

K -, W| + I " 21 < K (3.13) 

\W\ — W2\ 

for all wi,W2 € t(T) and for all w on the subarc of t(T) with smaller diameter. Thus circles and 
lines are 1-quasicircles. Quasicircles look very flat on all scales if K is close to 1, but for any K > 1 
they can contain a a dense set of spirals. See for example, Figure 8. 

If T satisfies (3.13) with K = 1 + S and small 5 and if {z k } C t(T) is locally evenly spaced 

then 



arg 



/ Z k - Zk-l 
\Zk+l — Zk 



(3.14) 
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for some constant C, depending on D. 

Theorem 3.11. There is a constant K > 1 so that ifT is a K -quasicircle with K = 1 + 5 < K 
and if {zk} are locally evenly spaced on T, then the geodesic algorithm finds a conformal map of 
H onto a region Q c bounded by a C (K)-quasicircle containing the data points {z k }. The constant 
C(K) can be chosen so that C(K) — > 1 as K — > 1. Moreover, given rj > 0, if the mesh size n({zk}) 
is sufficiently small then 

d H (T,dn c ) < V , 
where dn is the Hausdorff distance in the spherical metric. 

Proof. We may suppose that T satisfies (3.13) with K = 1 + 5 and 5 small. Note that oo G T. 
If {zk}\ are locally evenly spaced points on dfl, with /j, = max|zfc — Zk-i\ sufficiently small then 
(3.14) holds and D(oo, zi), D(z\, Z2), D(z n _i, z n ),D(z n , 00) is an C5^ -diamond chain, where the 
main axis of the cone D(oo, z\) is in the direction z\ — Z2 and the main axis of D(z n , 00) is in the 
direction z n — z n -\. Moreover D{oo, z\), D(zi, z 2 ), D(z n -i, z n ) satisfies the e-pacman condition 
if 

e > 

for some universal constant C. Now apply Theorem 3.4 to obtain jj C D(zj-i, Zj), j = 1, . . . , n — 1. 
By an argument similar to the proof of Theorem 3.10, we can also find a geodesic arc for C\(Uq _1 7j) 
from z n to 00 contained in D(z n , 00). Then the computed curve will be a C-fT-quasicircle. ■ 



As noted before, the boundary of the region computed with the geodesic algorithm, <9f2 c , is a 
C 1 curve. We end this section by proving that dQ c is slightly better than C 1 . If < a < 1, we say 
that a curve T belongs to C 1+a if arc length parameterization 7(5) of T satisfies 

W(S 1 )- 1 '(S 2 )\ <C\ 8l -8 2 \ a 

for some constant C < 00. 

We say that a conformal map / defined on a region Q belongs to C 1+a (fl), < a < 1 
provided / and /' extend to be continuous on Q and there is a constant C so that 

\f'(z 1 )-f'(z 2 )\<C\z 1 -z 2 r 



24 



for all zi, Z2 in U. 



Proposition 3.12. If the bounded Jordan region Q c is the image of the unit disc by the geodesic 
algorithm, then 

dn c g c 3 / 2 , 

and dQ c G" C 1+a for a > 1/2, unless Q c is a circle or a line. Moreover ip G C 3 / 2 (r2 c ) and 



Proof. To prove the first statement, it is enough to show that if 7 is an arc of a circle in EI which 
meets E orthogonally at (constructed by application of one of the maps f~ x as in Figure 2), 
then the curve a which is the image of [—1, 1] U 7 by the map S(z) = V z 2 — d 2 is C 5 (and no 
better class) in a neighborhood of S(0) = id. Indeed, subsequent maps in the composition ip~ x are 
conformal in M and hence preserve smoothness. For d > 0, the function 



2 



\ \1 + Vz 2 -c 2 /b / 



-dp = ld +^(z 2 - c 2 ) - ^(z 2 -c 2 )i+ 0((z 2 - c 2 ) 2 ) 



for some choice of b G E and c > is a conformal map of the upper half plane onto a region whose 

3 

complement contains the curve a. Clearly ip e near z = ±c, and so by a theorem of Kellogg 
(see [GM, page 62]), a G C*. The same theorem implies a is not in C a for a > | unless l/& = 0. 
This argument also shows that ip c G C 3 / 2 (r2). To prove (p~ l G C^(D), apply the same ideas above 
to the inverse maps. Alternative, this last fact can be proved by following the proof of Lemma 
II.4.4 in [GM]. ■ 



§4. Slits and Newton's method 

One complication of the slit and zipper algorithms is that the basic maps f a = g~ 1 cannot be 
written explicitly in terms of elementary maps, unlike the geodesic algorithm. Newton's method 
can be used to find the inverse of g a . 

Fix p, with < p < 1, and let f(z) = (z - p)P(z + 1 - pY~ p . Then /(H) =M\L where L is 
the line segment from to e inp p p (l - p) l ~ p . (Note that \ < \L\ < 1). Fix w G /(H). We wish to 
solve 

f(z)=w (4.1) 
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for z. Newton's method then takes an initial guess z and defines 

f(z n ) - w 



Zn+l Z-n 



f'(Zn) 



Near oo 



f(z) = z + l-2p + 0(- 



so a natural first guess for an approximation to the solution z to (4.1) would then be 



zq = w + 2p — 1. 



The next Theorem says that this initial guess zq will work for large w. 



Theorem 4.1. If < p < 1, set 



f( z ) = ( z -p)P( z + l-py-P, 



and suppose \w\ > (1 + \/5)/2. Then for z = w + 2p — 1 the n-th Newton iterate z n has relative 
error 



f(z n ) - w 



w 



< 



3 / 1 



2 V 12 



For example 



f(z 4 ) - w 



w 



< 10 



-17 



so that z 4 is virtually a formula for f~ 1 (w). In fact, in the slit or geodesic algorithm the points 
w with small \w\ correspond to points in the region near the corresponding vertex, so that most 
points will have large modulus. In practice, most points need only one or two iterations of Newton's 
method. The "approximate zero theorem" of Smale and Shub-Smale (see [SS]) can be also used to 
show that Newton's method will converge quadratically if |«;| is sufficiently large. Since we have 
an explicit formula for /, it is not surprising that we get a somewhat stronger result, in terms of 
\w\, in Theorem 4.1. 



Proof. Set F{z) = (f(z) - w)/w. We claim that 



\F(w + 2p-l)\ < 



p(l-p) 



and if 



\F(z)\ < 



p(l-p) 



(4.2) 



(4.3) 
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then 



\z\ 2 >p(l-p). 
To prove these claims, we study the auxiliary function 



H(C) = {1-(1-P)() P {1+P() 1 - P - 1, 



which has derivative 



p(p- i)C 



[i + (p-i)C] 1 - p [i+pC] p ' 

Bounding the denominator from below and integrating we obtain the estimate 

p(l-p) |C| 2 

|F(C)I * 2 i-icr 

Note first that F(w + 2p-l) = H(l/w). So that by (4.5) 

+V5 



2 V1 + V57 1- T 



proving (4.2). 

Suppose now that (4.3) holds and \z\ 2 < p(l —p). Then 



(z - p)P(z + 1 - pf-P < p(l - p) 1 



This implies 



it; 



8 i i i i 

< ^[P 5 (l - f) 5 +P] P [P 5 (1 - p) 2 + 1 - 



y[p p (l-p) 1 - p (l + 2p'(l-p)')]' 



8\/2 y^+l 
" 7 < 2 ' 

contradicting our assumption \w\ > (\/5 + l)/2, and proving that (4.3) implies 
Next suppose that (4.3) holds and set 

' z + l-p xP 



z = z - ^V^r = -2 + ' 



z — p 



F'(z) \ z 

Then after some manipulations we obtain the magic formula 

Fiz)= H ( F ^ 



-(z + l-p) 
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By (4.3) and (4.4), 



F(z) 



yjp{\-p) 1 

2 - 4' 



and so by (4.6), (4.5), (4.4) and (4.3) 



By induction and (4.2) 



proving Theorem 4.1. 



W2)l < 



p(l -p) 



2(S) 



<l\F(z)f< P -^ 



§|F<*.)| < (|f (2 „)|) < (1 N 



The region of possible w where quadratic convergence is obtained can be enlarged with more 
involved estimates, but Newton's method applied directly to / with this initial value will not always 
converge. Indeed, the Newton interate z — F(z)/F'(z) has repelling fixed points at p and p — 1, 
and a pole at 0. In order to successfully apply the algorithm to a wide variety of curves we need to 
find a reliable routine for finding the inverse. 

In the implementation of the slit and zipper algorithms we consider four regions based on the 
length \L\ = p p (l — p) 1_p of the segment L and the imaginary part of the tip w t i p = /(0). 



^oo = {w : \w\ > -\L\} 



^ti P = {w :\w - w ti p | < -Imwtip} 



Q p = {w : < arg w < irp} 



{w : up < arg w < ir} 
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If w € f^oo then we use Newton's method as described in Theorem 4.1. However, we improve 
the initial guess by taking more terms in the expansion at oo: 

L0 -i | P(l ~P) . (1 -2p)p(l -p) 
zo = u; + 2p-l + ^^ + — , 

and we rewrite the function to iterate on z/w instead of z to improve numerical accuracy. If 
w $l f^oo but w G fitip then we first open up the region by applying k(w) = y/w — w t \ p . Then k o f 
extends to be analytic and one-to-one in a neighborhood of 0. So we use Newton's method to solve 
k o f(z) = k(w) for z. The remaining w are in the sectors between R and L. If w $ ft^ U Q, t - ip 
but < argtt; < np, then we apply the preliminary map k p (z) = instead of k and use Newton's 
method again. For the remaining points we use the preliminary map ki— p (z) = z 1 ^ . We leave 
the proof of the analog of Theorem 4.1 in the remaining three regions to the interested reader. 
While we cannot prove convergence of Newton's method in every case, extensive numerical testing 
indicates that we have chosen the proper regions. 



§5. Estimates for conformal maps onto nearby domains 

We begin this section with a discussion of the following question. Consider two simply con- 
nected planar domains f2j with £ Qj and conformal maps ifj : Q,j — > D fixing 0, suitably 
normalized (for instance positive derivative at 0). If Q,\ and 2 are "close," what can be said about 
\(fi —(f2\ on J7i (1^2, or about [ip^ 1 — i^ 1 ) on D? The article [W] gives an overview and numerous 
results in this direction. How should "closeness" of the two domains be measured? Simple examples 
show that the Hausdorff distance in the Euclidean or spherical metric between the boundaries does 
not give uniform estimates for either \\ipi — c^lloo or — ^ 1 || 00 . 




Figure 11. Small Hausdorff distance 
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For example in Figure 11, f2i contains a disc of radius 1 — 8 where 8 is small and hence for 
VL 2 = B, dff(fii,02) < 5, but |<^i(-Zi) — vpi^)! is large and |<^i(^2) — ^1(^3) I is small so that 
neither ||¥>i(z) — z\\oo nor — z\\oo is small. 

Mainly for ease of notation, we will assume throughout this section that the Qj are Jordan 
domains, and denote 7j : 9B — ► dQ.j an orientation preserving parametrization. Even the more 
refined distance 

inf H71 - 72 o alloc, 

a 

where the infimum is over all homeomorphisms a of <9B, does not control — ( p 2 1 \\ 00 or \ \ipi — 

<^2||oo- For example, let VL 2 be a small rotation of the region in Figure 11. What is needed is 
some control on the "roughness" of the boundary. Following [W], for a simply connected domain 
O we define 

v($) = Vn(8) = supdiam T(C), 

c 

where the supremum is over all crosscuts of £1 with diam C < 8, and where T(C) is the component 
of Q \ C that does not contain 0. Notice that rf(8) — ► as 8 — > is equivalent to saying that dQ is 
locally connected, and the condition rj(8) < K8 for some constant K is equivalent to saying that 
VL is a John-domain (e.g. [P], Chapter 5). It is not difficult to control the modulus of continuity 
of ip -1 : B — ► $7 in term of rj, see [W], Theorem I. This can be used to estimate \ \ipi 1 — ^^lloo i n 
terms of the Hausdorff distance between the boundaries, for example. 

Theorem 5.1 (Warschawski[W], Theorem VI). If Qi and £l 2 are John-domains, i]j{8) < n8 
for j = 1, 2, and if du(dfli, dfl 2 ) < e, then 

llvT'-^lloo <Ce<* 

with a = a(K) and C = C(k, dist(0, 0Qi U dQ 2 )). 

In fact, Warschawski proves that every a < 2/(tt 2 k?) will work (with C = C(a)). Using the 
Holder continuity of quasiconformal maps, his proof can easily be modified to give the following 
better estimate if Oi and 2 are i^-quasidiscs with K near 1. A K-quasidisc is a Jordan region 
bounded by a i^-quasicircle. 

Corollary 5.2. IfQx and £l 2 are K-quasidiscs, and if d/j(<9f2i, Sf^) < e, then 

Wvl 1 ~ V^Woo <Ce a 
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with a = a(K) — > 1 as K — > 1. 

As for estimates of ||^>i — 992II00, Warschawski shows [W, Theorem VII] that 

2 

sup — 9^2 1 < Ce 1 / 2 log - 
fii e 

if J7i C f2 2 , and if f2i is a John-domain, with C depending on k and on dist(0, <9f^i UdCl 2 ). However, 
his result does not apply without the assumption of inclusion Hi C ^2. To treat the general case 
the trick of controlling \<pi — <p 2 \ by passing to the conformal map ip of the component O of f^i n £l 2 
containing (which now is included in Qj) does not seem to work, as the geometry of ft can not 
be controlled. Nevertheless, for the case of disc-chain domains, the above estimate can be proved, 
even without any further assumption on the geometry on the circle chain: 

Theorem 5.3. Let D 1: D 2 , .-.,D n be a closed e-disc-chain surrounding 0. Suppose dQj C L>kDk 
for j = 1,2, and let ipj : Qj — > D be conformal maps with ipi (0) = 992 (0) = and (f\ (p) = ip 2 (p) for 
a point p 6 <9f2i PI <9S1 2 - Then 

SUp \ifi 

where C depends on dist(0, UkDk) only. 

In case we have control on the geometry of the domains, we have the following counterpart to 
Corollary 5.2. 

Theorem 5.4. If and £l 2 are K-quasidiscs, if du{d£li, d£l 2 ) < e, and if tpi(pi) = ^2(^2) for a 
pair of points pj £ dQj with \pi — p 2 \ < e, then 

sup \<pi(w) - (f 2 (w)\ < Ce a 
wen 

with a = a(K) — > 1 as K — > 1, where Q, is the component of Q\ n ^2 containing 0. 

The proofs of both theorems rely on the following harmonic measure estimate, which is an 
immediate consequence of a theorem of Marchenko [M] (see [W, Section 3] , for the statement and a 
proof). To keep this paper self-contained, we include a simple proof, shown to us by John Garnett, 
for which we thank him. 

Lemma 5.5. Let < B < it, < e < 1/2 and set D = D \ {re u : -6 < t < 9, 1 - e < r < 1}, 
A = dD\ dB. Then 

w(0,A,D) < - + Celog- 
7r e 
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for some universal constant C. 



Proof. Set to(z) = u(z,A,D) for z G D. By the mean value property, it is enough to show that 

<*) < c^- e 

for z = (1 — e)e rf and + e < i < it. To this end, set I = {e lT : —6 < r < 9} and consider the 
circular arc {(" : uj((,,I, D) = |}. If e < eg for some universal £q (for e > eo there is nothing to 
prove), then A is disjoint from this arc and it follows that oj((,I, B) > | on A. The maximum 
principle implies w(C) < 3u;(C, i", B) on D. Now the desired inequality follows from 

u((l-e)e it ,I,0) = — [ - 1 lo dr < Ce / . 1 .„ dT < C— — n . ■ 

U ^ ^ 2vr7_ e |(l- e ) e .*- e *r|2 - J_ e (t- T )2 t _ e 

Proof of Theorem 5.3. We may assume that <Pj(j?) = 1. We will first assume that p is one of the 
points Dfefl-Dfe+i. Denote the largest simply connected domain C C containing whose boundary 
is contained in Ufc-Dfc (thus O is the union of UkDk and the bounded component of C \ UkDk), and 
(p the conformal map from to D with <^(0) = and <p(p) = 1. First, let z £ dCli n 90. Denote 
-B respectively i?i the arc of d£l (dQi) from p to z. By the Beurling projection theorem (or the 
distortion theorem), every <p(Dj) has diameter < C\fe. Therefore ip(B\) is an arc in B, with same 
endpoints as <p(B), that is contained in S = {re lt : 1 — < r < 1, —C^fe < t < dxgip(z) + C^fe\. 
Denote ^4 = By Lemma 5.5, 

w(0,Bi,ni) < w(0,Bi,fi\Si) < w(0,A,D\A) < — argp(z) + 2(7^ + Cy^log -7= 

2-k yje 

and we obtain 

wgip^z) = 27Ttx;(0,Bi,ni) < arg(^(z) + Ce 1/2 log-. 

The same argument, applied to the other arc from p to 2, gives the opposite inequality, and together 
it follows that 

\<p(z)-<pi(z)\ <Ce^ 2 \og^. 
Now let z G be arbitrary. If z' is a point of dCli n 90 in the same disc Dj as z, then we 

have 

Hz) - <fi(z)\ < \<p(z) - <p(z')\ + \<p(z') - + \<Pi(z) - Mz')\ < 2CV~e + Ce 1 ' 2 log \. 
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The maximum principle yields \ip — ipi\ < Ce 1 / 2 log \ on f^. The same argument applies to \ip — ip 2 \, 
and the theorem follows from the triangle inequality. 

If p 6 dQiHd^2 is arbitrary, let p' be one of the points DkCiDk+i in the same disc Dj as p. Then 
the above estimate, applied to a rotation of ipi, ip2 and p' gives \<P2(p')/<Pi(pf)<Pi — V2I < Ce 1 / 2 log 2 
and the theorem follows from |y>j(p) — <Pj(p')\ ^ C\/e. ■ 

The following lemma is another easy consequence of the aforementioned theorem of Marchenko [M] 
([W], Section 3). 

Lemma 5.6. Let H C B be a K-quasidisc with £ H such that dH C {1 — e < \z\ < 1}, and 7et 
h be a conformal map from B to H with h(0) = and \h(p) — p\ < e for some p 6 <9B. Then 

\h(z) -z\< Celog^, 

where C depends on K only. 

Proof. We may assume p = 1. Let 2 = e lT and consider the arc A = {/i(e lt ) : < t < r} C dif 
of harmonic measure t/2tt. For suitable C = C(K) we have that D = B \ {re lt : — Ce < t < 
argh(z) + Ce, 1 — e < r < 1} contains A. By the maximum principle and Lemma 5.5, 

r/2vr = u(0, A,H) < u(0, dD n B, D) < aigh(z)/2n + Celog i. 

Applying the same reasoning to <9if \ A, the lemma follows for all z 6 <9B and thus for all zeD.I 

Note that the conclusion of Lemma 5.6 is true if instead of assuming if is a i^-quasidisc, we 
only assume arg z is increasing on dH. 

Proof of Theorem 5.4. Because fii and Q2 are K-quasidiscs, ipi and ip 2 have if 2 -quasiconformal 
extensions to C (see [L], Chapter 1.6). In particular, they are Holder continuous with exponent 
1/K 2 (see [A]), and it follows that with a = 1/K 2 and r = 1 - Ce a , we have ^({^l ^ r l) c ^2- 
In particular, h(z) = ip 2 ((Pi 1 (rz)) is a conformal map from B onto a iT 4 -quasidisc H C B, and by 
the Holder continuity of (p 2 and ip^ 1 we have dH C {1 — Ce" 3 < \z\ < 1}. Now Lemma 5.6 yields 
\h(z) - z\< Ce 13 , for any (3 < a 3 and C = C((3). For w £ Q C £li n ft 2 , let z = <pi(w), then 

biH - ¥>2(w)i = \z- ^(^r 1 ^))! < \ z - ^farV*))! + 1^2 (^r 1 ^)) - ^(^r 1 ^))! < c e/3 > 
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where again we have used the Holder continuity of if2 and <p 1 . The Theorem follows. 



§6. Convergence of the Mapping Functions 

We will now combine the results of Sections 2 and 3 with the estimates of the previous section, 
to obtain quantitative estimates on the convergence of the geodesic algorithm. Throughout this 
section, Q will denote a given simply connected domain containing 0, bounded by a Jordan curve 
dQ, z , z n are consecutive points on d£l, Q c is the domain and ip c : O c — > D the map computed 
by the geodesic algorithm, and ip : Q — > ID is a conformal map, normalized so that ip c (0) = <p(0) = 
and <Pc{po) = <f(Po) for some po € <90 n <9fi c . 

Combining Theorems 2.2 and 5.3 and Propositions 2.5 and 3.12 we obtain at once: 

Theorem 6.1. If <9$7 is contained in a closed e-disc-chain Uj=o an< ^ ^ z i = ^ then 
<9f2 c is a smooth (C^) piecewise analytic Jordan curve contained in \J™ =0 Dj U Zj, the map ip c 
extends to be conformal on Q U £l c and 



Now assume that dfl is a i<C-quasicircle with K < K and assume approximate equal spacing 



where d (essentially the Minkowski-dimension) is close to 1 when K is close to 1. Combining 
Theorem 3.11 with Corollary 5.2 and Theorem 5.4, we have: 

Theorem 6.2. Suppose <9f2 is a K -quasicircle with K < Kq. The Hausdorff distance between 
d£l and dfl c is bounded by C(K)e, where C(K) tends to as K tends to 1 and n to infinity. 
Furthermore, 



sup \ip(w) — <p c (w)\ < Ce 1 / 2 log -. 
wen ' e 



of the Zj, say, |e < \zj + \ — Zj\ < 2e. Then 



C C 
— < n < -r 
e ~ ~ e d 



(6.1) 



<P 1 - 



-l 




< Ce 



and 



sup \<p(w) — ip c (w)\ < Ce 1 



ce 
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with a = a(K) — > 1 as K — > 1, where JIq is the component of n fi c containing 0. 



The best possible exponent in (6.1) in terms of the standard definition of if(<90), which 
slightly differs from our geometric definition, is given by Smirnov's (unpublished) proof of Astala's 
conjecture, 

This allows us to easily convert estimates given in terms of e, as in Theorem 6.2, into estimates 
involving n. 

Finally, assume that dQ. is a smooth closed Jordan curve. Then £1 is a f^T-quasicircle and a 
John domain by the uniform continuity of the derivative of the arc length parameterization of dQ. 
The quasiconformal norm K(dQ) and the John constant depend on the global geometry, as does 
the e-pacman condition when there are not very many data points. As the example in Figure 11 
shows, even an infinitely differentiable boundary can have a large quasiconformal constant and a 
large John constant. However, the e-pacman condition becomes a local condition if the mesh size 
^{{zk}) = maxfc \zk+i — Zk\ of the data points is sufficiently small. The radii of the balls in the 
definition of the e-pacman condition 

R k = C 1 lZk+1 ~ Zkl (6.2) 

increase as e decreases, but can be chosen small for a fixed e if the mesh size fi is small. To 
apply the geodesic algorithm we suppose that the data points have small mesh size and, as in the 
proof of Theorem 3.10, \(zq — z n )/(z n -\ — z n )\ is sufficiently large so that the e diamond chain 
D(zo, zi), . . . , D(z n -i, z n ) satisfies the e-pacman condition and 

n 

,Zk+l) 

k=0 

where D(z n , z n+ {) = D(z n ,z ) is an e-diamond. This can be accomplished for smooth curves by 
taking data points z , . . . , z n , zq with small mesh size and discarding the last few z n - ni , ■ ■ ■ ,z n 
where m is an integer depending on e and on dQ. The remaining subset still has small mesh size 
(albeit larger). This process of removing the last few data points is necessary to apply the proof of 
Theorem 3.10, but in practice it is omitted. We view it only as a defect in the method of proof. 

If dQ G C 1 and if tp is a conformal map of Q onto B then arg((^ -1 )' is continuous. Indeed, 
it gives the direction of the unit tangent vector. However there are examples of C 1 boundaries 
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where tp' and (v? -1 )' are not in continuous. In fact it is possible for both to be unbounded. If we 
make the slightly stronger assumption that dfl G C 1+a for some < a < 1, then ip G and 
p~ x € by Kellogg's theorem (see [GM, page 62]). In particular the derivatives are bounded 

above and below on f2 and D, respectively. Because of Proposition 3.12, we will consider the case 
1 + a = 3/2. Similar results are true for 1 + a < 3/2. 

Theorem 6.3. Suppose <9fi is a closed Jordan curve in C 3 / 2 and (p is a conformal map of onto 
D. Suppose zo, Zi, —z n , zq are data points on with mesh size fx = max \zj — Zj + i\. Then there 
is a constant C\ depending on the geometry of d£l, so that the Hausdorff distance between 90 and 
dQ, c satishes 

d H (dn,dn c ) < d^' 2 (6.3) 

and the conformal map ip c satishes 

\\<p- 1 -Vc 1 \\°o<Cn. p (6.4) 

and 

sup \p(z)-p c (z)\<C^, (6.5) 

for every p < 3/2. 



For example if n data points are approximately evenly spaced on dtl, so that fi = C/n then 
the error estimates are of the form C/n 3 / 2 in (6.3) and C/n p for p < 3/2 in (6.4) and (6.5). While 
Theorem 6.3 gives simple estimates in terms of the mesh size or or the number of data points, 
smaller error estimates can be obtained with fewer data points if the data points are distributed 
so that there are fewer on subarcs where d£l is flat and more where the boundary bends or where 
it folds back on itself. In other words, construct diamond chains with angles satisfying the 
£fc-pacman condition centered at Zk for each k. The errors will then be given by 

max I E k \z k - z k+ i 




Proof. It is not hard to see from (6.2) that dQ. satisfies the e-pacman condition with 
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for C sufficiently large. By the proof of Theorem 3.10, <9f2 c is contained in the union of the 
diamonds. The diamonds D(zk, Zfc+i) nave angle C/i 1 / 2 and width bounded by C/i and therefore 
(6.3) holds. 

Let tp be a conformal map of ID onto the complement of S7, C* \ fl. Then by Kellogg's Theorem 
as mentioned above, ip G C 3 / 2 . In particular, |^/| is bounded above and below on 1/2 < \z\ < 1. 
By the Koebe distortion theorem there are constants Ci, C2 so that 

Ci(l - |z|) < dist(^(«),SS2) < C 2 (l - \z\), 

for all 2 with 1/2 < |z| < 1. Thus we can choose r = 1 — C3/X 3 / 2 so that the image of the circle of 
radius r, I r = ip({\z\ = r}), does not intersect the diamond chain and da (In dfl) ~ [i 3 / 2 . Then the 
bounded component of the complement of I r is a Jordan region U r containing 0, and bounded by 
I r e C 3 / 2 , with C 3 ' 2 norm dependent only on dCl, and the bounds on \ip'\. 

Let <j be a conformal map of U r onto B. Inequality (6.4) now follows from [W, Theorem VIII] 
by comparing the conformal maps ip~ x and if' 1 to the conformal map <7 _1 where a : U r — > D and 
where all three (inverse) conformal maps are normalized to have positive derivative at and map 
to the same point in Q. 

To see (6.5), note that 

a(dn U dn c ) C{z:l-\z\< c/i 3/2 }. 

Moreover, because dfl U d£l c is contained in the diamond chain, and because both a G C 3 ! 2 and 
<7 _1 G C 3 / 2 , argcr(C) is increasing along dQ, for ^ sufficiently small. By the remark after the proof 
of Lemma 5.6, 

|w(0,7,<r(n)) -w(0,7*,B)| < C/x 3 / 2 ^^ 

for every subarc 7 of <t(9$7), where 7* denotes the radial projection of 7 onto dO. The same 
statements are true for dVi c . Then (6.5) follows because the harmonic measure of the subarc 7 P of 
dfl from p to p is given by 

and a similar statement is true for ^> c . ■ 

The constant C in Theorem 6.3 depends on the quasiconformality constant K(dQ), p, diam(fi), 
dist(0, dfl), and on 

M= sup (h//|,l/h//|), 

1/2<|z|<1 
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where ip is a conformal map of the complement of O to D. If / r = ^({M = r }) is replaced by 
a C 3 / 2 curve which is constructed geometrically instead of using the conformal map ip, then the 
constant C can be taken to depend only on the geometry of the region Q. 

§7. Some Numerical Results 

An in depth comparision of the algorithms in this article with other methods of conformal 
mapping and convergence rates will be written separately. To give the reader a sense of the speed 
and accuracy of computations, if 10,000 data points are given, it takes about 20 seconds with the 
geodesic algorithm to compute the mapping functions on an 3.2 GHz Pentium IV computer. Since 
all of the basic maps are given explicitly in terms of elementary maps, the speed depends only 
on the number of points, not the shape of the region or the distribution of the data points. The 
accuracy can be measured if the true conformal map is known. For example 

1 + yrzy 

where r < 1 maps the unit disc into an inverted ellipse. See Figure 12. The region was chosen be- 
cause it almost pinches off at 0, and because the stretching/compression given by max |/'|/min |/'| 
is big for r near 1. Higher resolution images can be obtained from: 

|http: / /www.math.washington.edu / ^--marshall / preprints /zipper .pdf 




Figure 12. Inverted ellipse with r = .95. 



We chose r = .95 and used as data points the image by / of 10,000 equally spaced points on the unit 
circle, and compared the corresponding points on the unit circle computed by the geodesic algorithm 
with 10,000 equally spaced points. The errors were less than 1.8 • 10 -6 . The same procedure using 
the zipper algorithm took 84 seconds, and had errors less than 9.2 • 10 -8 . When the number of data 
points was increased to 100,000, the time to run the geodesic algorithm increased to 25 minutes 
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with errors less than 2 • 10 -8 . In this example, the the difference between successive boundary data 
points ranged from .025 to 3 • 10 -6 so that perhaps a better distribution of data points would have 
given smaller errors. 

Figure 13 shows the conformal map of a Carleson grid on the disc to both the interior and 
exterior of the island Tenerife (Canary Islands). The center of the interior is the volcano Teide. 
It also shows both the original data for the coastline, connected with straight line segments, and 
the boundary curve connecting the data points using the zipper algorithm. At this resolution, it is 
not possible to see the difference between these curves. The zipper algorithm was applied to 6, 168 
data points and took 36 seconds. The image of 24, 673 points on the unit circle took 48 seconds 
and all of these points were within 9 • 10 -5 of the polygon formed by connecting the 6, 168 data 
points. The points on the circle corresponding to the 6, 168 vertices were mapped to points within 
10 -10 of the verticies. This error is due to the tolerence set for Newton's method, round-off error, 
and the compression/expansion of harmonic measure. The image of 8, 160 verticies in the Carleson 
grid took 25 seconds to be mapped to the interior and 25 seconds to the exterior. 



Figure 13. Tenerife. 

Higher resolution images can be obtained from: 

http: / /www. math.washington.edu/~marshall / preprints /zipper .pdfj 

The first objection one might have in applying these algorithms with a large number of data 
points is that compositions of even very simple analytic maps can be quite chaotic. Indeed this is 
the subject of the field complex dynamics. We could redefine the basic maps f a by composing with 
a linear fractional transformation of the upper half plane so that the composed map is asymptotic 
to z as z — > oo. This will not affect the computed curve in these algorithms since the next basic 
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map begins with a linear fractional transformation (albeit altered). However, if we formulate the 
basic maps in this way, then because the maps are nearly linear near oo, the numerical errors will 
accumulate only linearly 

Osculation methods also approximate a conformal map by repeated composition of simple 
maps. See Henrici [H] for a discussion of osculation methods and uniform convergence on compact 
sets. The algorithms of the present article follow the boundary of a given region much more closely 
than, for instance, the Koebe algorithm and give uniform convergence rather than just uniform 
on compacta. It is possible to use the techniques of this paper to prove the geodesic algorithm is 
an osculation method for smooth curves, and therefore by the results in [H] converge uniformly 
on compact subsets. However, prior to this article even a proof that these methods satisfied the 
osculation family conditions was not known. 

Recently Banjai and Trefethen [BT] adapted multigrid techniques to the Schwarz-Christoffel 
algorithm and successfully computed the conformal map to a region bounded by a polygon with 
about 10 5 edges. They used a 12 fold symmetry in the region to immediately reduce the parameter 
problem to size 10 4 . Any other conformal mapping technique can also use symmetry and obtain a 
12 fold reduction in the number of data points required, however their work does show at least that 
Schwarz-Christoffel is possible with 10 4 vertices, though convergence of the algorithm to solve the 
parameter problem is not always assured. The zipper algorithm is competitive in speed and accuracy 
for such regions. The geodesic algorithm is almost as good, and has the advantage that it is very 
easy to code and convergence can be proved. It would be interesting to try to prove convergence of 
the technique used in [BT] to find the prevertices, for polygons which are -fT-quasicircles in terms 
of K. It would be interesting as well to apply multigrid techniques to the zipper algorithm. 

One additional observation worth repeating in this context is that the geodesic and zipper 
algorithms always compute a conformal map of EI to a region bounded by a Jordan curve passing 
through the data points, even if the disc-chain or pacman conditions are not met. The image region 
can be found by evaluating the function at a large number of points on the real line. By Proposition 
2.5 and Corollary 3.9, if the data points {zj} satisfy the hypotheses of Theorem 2.2 or Theorem 
3.4, then ip can be analytically extended to be a conformal map of the original region f2 to a region 
very close to B. To do so requires careful consideration of the appropriate branch of \fz at each 
stage of the composition. 

Theorem 2.2 and Theorem 3.4 and their proofs suggest how to select points on the boundary 
of a region to give good accuracy for the mapping functions. Roughly speaking, points need to 
be chosen closer together where the region comes close to folding back on itself. See Figure 12 
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for example. Greater accuracy can be obtained by placing more points on the boundary near the 
center and fewer on the big lobes. See also the remarks after the statement of Theorem 6.3 in this 
regard. In practice, the zipper map works well if points are distributed so that 

B(z k ,5\z k+1 -z k \) ndtt (7.1) 

is connected. 

When the boundary of the given region is not smooth, then one of the processes described in 
section 2 should be used to generate the boundary data, if the geodesic algorithm is to be used. 
For example, if nothing is known about the boundary except for a list of data points, then we 
preprocess the data by taking data points along the line segments between the original data points, 
so that these new points correspond to points of tangency of disjoint circles centered on the line 
segments, including circles centered at the original data points. Note that the original boundary 
points are not among these new data points. The geodesic algorithm then finds a conformal map to 
a region with the new data points on the boundary. The boundary of the new region will be close 
to the polygonal curve through the original data points, but will not pass through the original data 
points. This boundary is "rounded" near the original data points. Indeed it is a smooth curve. 

When the boundary of the desired region is less smooth, for example with "corners" , then the 
zipper or slit algorithms should be used. In this case additional points are placed along the line 
segments between the data points, with at least 5 points per edge and satisfying (7.1). In practice, 
at least 500 points are chosen on the boundary so that the image of the circle will be close to the 
polygonal line through the data points. Since two data points are pulled down to the real line with 
each basic map in the zipper algorithm, the original data points should occur at even numbered 
indices in the resulting data set (the first data point is called Zq). Then the computed boundary 
Q c will have corners at each of the original data points, with angles very close to the angles of the 
polygon through the original data points. 

A version of the zipper algorithm can be obtained from [MD]. The conformal mapping programs 
are written in Fortan. Also included is a graphics program, written in C with X-ll graphics by 
Mike Stark, for the display of the conformal maps. There are also several demo programs applying 
the algorithm to problems in elementary fluid flow, extremal length and the hyperbolic geometry. 
Extensive testing of the geodesic algorithm [MM] and an early version zipper algorithm was done 
in the 1980's with Jim Morrow. In particular that experimentation suggested the initial function 
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<po in the zipper algorithm which maps the complement of a circular arc through zq, Z\, and z 2 
onto H. 



Appendix. J0rgensen's Theorem. 

Since J0rgensen's theorem is a key component of the proof of the convergence of the geodesic 
algorithm, we include a short self-contained proof. It says that discs are strictly convex in the 
hyperbolic geometry of a simply connected domain £1 (unless dtt is contained in the boundary of 
the disk). 

Theorem A.l (J0rgensen [J]). Suppose Q is a simply connected domain. If A is an open disc 
contained in f2 and if 7 is a hyperbolic geodesic in Q, then 7 n A is connected and if non-empty, it 
is not tangent to <9A in Q. 



Proof. See [P, page 91-93]. Applying a linear fractional transformation, we may suppose that the 
upper half plane H C tt. Suppose x € R and suppose that / is a conformal map of D onto Q such 
that /(0) = x and /'(0) > 0. Then 

is a bounded harmonic function on D which is greater than or equal to by the maximum principle. 
Thus ImjQ^ > on ( — 1,1) and hence Im/(2;) < on the diameter (—1,1). The condition 
/'(0) > means that the geodesic /((— 1, 1)) is tangent to R at x. Thus if 7 is a geodesic which 
intersects HI and contains the point x, then it cannot be tangent to R at x. Two circles which are 
orthogonal to R can meet in EI in at most one point, and hence hyperbolic geodesies in simply 
connected domains (images of orthogonal circles) meet in at most one point. Thus 7 cannot reenter 
EI after leaving it at x because it is separated from R by the geodesic /((— 1, 1)). The Theorem 
follows. ■ 

In Section 2, we commented that a constructive proof of the Riemann mapping theorem fol- 
lowed from the proof of Theorem 2.2. The application of J0rgensen's theorem in the proof of 
Theorem 2.2 is only to domains for which the Riemann map has been explicitly constructed. 
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